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We study analytically defect polariton states in Bragg multiple-quantum-well structures and defect 
?H . induced changes in transmission and reflection spectra. Defect layers can differ from the host layers in 

three ways: exciton-light coupling strength, exciton resonance frequency, and inter-well spacing. We 
show that a single defect leads to two local polariton modes in the photonic bandgap. These modes 
cause peculiarities in reflection and transmission spectra. Each type of defect can be reproduced 
experimentally, and we show that each of these plays a distinct role in the optical properties of 
the system. For some defects, we predict a narrow transmission window in the forbidden gap at 
the frequency set by parameters of the defect. We obtain analytical expressions for corresponding 
local frequencies as well as for reflection and transmission coefficients. We show that the presence 
of the defects leads to resonant tunneling of the electromagnetic waves via local polariton modes 
^ ■ accompanied by resonant enhancement of the field inside the sample, even when a realistic absorption 

is taken into account. On the basis of the results obtained, we make recommendations regarding 
the experimental observation of the effects studied in readily available samples. 
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I. INTRODUCTION 



Optical properties of multiple quantum wells (MQW) have attracted a great deal of interest recently.0"! Unlike 
' other types of superlattices, excitons in MQW are confined in the planes of the respective wells, which are separated 
by relatively thick barriers. Therefore, the only coupling between different wells is provided by the radiative optical 
' field. The coupling results in MQW polaritons - coherently coupled quasi-stationary excitations of quantum well 
(QW) excitons and transverse electromagnetic field. The spectrum of short MQW structures consists of a number of 
quasi-stationary (r.adiative) modes with finite life-times. This spectrum is conveniently described in terms of super- or 
sub-radiant modesEmti When the number of wells in the structure grows, the life-time of the former decreases, and the 
life-time of the latter increases. In longer MQW structures, however, this approach becomes misleading, as discussed 
in Ref. ^, and a more appropriate description is obtained in terms of stationary modes of an infinite periodic structure. 
The spectrum of MQW polaritons in this .case consists of two branches separated by a gap with a width proportional 
to the exciton-light coupling constant Totl The length at which this description becomes preferable, usually depends 
upon a problem at hand. For instance, even though the discrete structure of the sub-radiant modes of 100 long 
I ■ MQW in Refs. MM could be resolved, the description of the gap region in terms of the super-radiant mode leads to an 
apparent contradiction with the absence of significant luminescence in this region. If one is interested in properties 
of the gap region, a "long system" usually means that it is longer than the penetration length of radiation into the 
sample. The latter length depends upon the frequency, therefore the system can be long enough for frequencies close 
to the center of the gap, and still "short" for frequencies in the vicinity of the band edges. Systems similar to those 
studied in Ref. ^ can be considered sufficiently long for 90% of the gap region when the number of wells becomes close 
to 200. nnnn 

In a number of papers3l§LK3 it was shown that the width of the polariton bandgap can be significantly increased 
& \ by tuning the interwell spacing, a, to the Bragg condition, a = Ao/2, where Ao is the wavelength of the light at 
the exciton frequency Qq. Under this condition, boundaries of two adjacent gaps become degenerate, and one wider 
gap with the width proportional to VT is formed. Detuning of the lattice constant from_the exact Bragg condition 
removes the degeneracy and gives rise to a conduction band in the center of the Bragg gapJj A well-pronounced Bragg 
polariton gap was observed in recent experimentscl with InGaAs/GaAs Bragg structures with the number of wells 
up to 100. These experiments convincingly demonstrate that despite homogeneous and inhomogeneous broadening, 
the coherent exciton-photon coupling in long MQW is experimentally feasible. Polariton effects arising as a result of 
this coupling open up new opportunities for manipulating optical properties of quantum heterostructures. 

One such opportunity is associated with introducing defects in MQW structures. These defects can be either 
QW's of different compositions replacing one or several "host" wells, or locally altered spacing between elements of 
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the structure. It is well known in the physics of regular crystals (see, for instance, Ref. [10[) that local violations of 
otherwise periodic structures can lead to the appearance of local modes with frequencies within spectral gaps of host 
structures. This idea was first applied to MQW in Ref. [□], where it was shown that different defects can indeed give 
rise to local exciton-polariton modes in infinite MQW. Unlike regular local modes in 3-d periodic structures, these 
modes are localized only in the growth direction of the MQW, while they can propagate along the planes of the wells. 
Therefore, one should clearly distinguish these defect polariton modes fpara_well known interface modes in layered 
systems or non-radiative two-dimensional polariton modes in ideal MQW.II3~li-a The latter exist only with the in-plane 
wave numbers, fen, exceeding certain critical values, while the local mode in a defect MQW structure exists at fen = 
and can be excited even at normal incidence. 

In this paper we present results of detailed studies of local polariton modes (LPM's) produced by four different 
types of individual defects in Bragg MQW structures. The peculiar structure of Bragg MQW's results in a wider than 
usual polariton gap, which is actually formed by two gaps with degenerate boundaries. This property has a profound 
effect on the properties of local modes, leading, for instance, to emergence of two local states from a single defect. 
The similar effects of doubling the number of local modes was also predicted in the case of braggoritons - excitations 
arising inside photonic band gaps of periodic structures made of a dispersive material, if the resonance frequency of 
the dielectric functions happens to belong to the gap.ta 

We neglect in-plane disorder in individual QW's, and assume that apart from a deliberately introduced single defect, 
the structure remains ideally periodic. We consider LPM's with zero in-plane wave vector only. Such modes can be 
excited by light incident in the growth direction of the structure, and can result in the resonance transmission of light 
with gap frequencies. This effect is studied both analytically and numerically with homogeneous broadening taken 
into account phenomenologically. Equations describing dynamics of the light-exciton interaction in this situation are 
essentially equivalent to a model of the oMpdimensional chain of dipoles used in the series of our previous works whepa 
LPM's in polar crystals were discussed. Similar equations also appear in the theory of atomic optical lattices.^ 
The essential difference between results presented in this paper and previous studies stems from the peculiarities of the 
Bragg arrangement. Using Greens' function and transfer matrix formalisms, we study both eigenfrequencies of LPM's 
for different types of defects and transmission properties of the defect structures. Taking into account homogeneous 
broadening and using parameters of the system studied experimentally in Refs. ^,^|, we predict which defects will 
produce the most significant changes in transmission and reflection properties of realistic MQW structures. 



In order to describe optical properties of QW's one has to take into account the coupling between retarded electro- 
magnetic waves and excitons. This is-usually done with the use of the non-local susceptibility determined by energies 
and wave functions of a QW exciton.Elo The treatment of the exciton subsystem can be significantly simplified if the 
interwell spacing is much larger than the size of a well itself. In InGaAs/GaAs MQW structures studied in Ref. ||, 
on which we base our numerical examples, the width of the QW layer amounts only to about 10% of the period of the 
structure. In this case, one can neglect the overlap of the exciton wave functions from neighboring wells and assume 
that an interaction between well excitons occurs only due to coupling to the light. It is also important that the width 
of the wells is also considerably smaller than the exciton's Bohr radius, and, therefore, one can neglect spatial extent 
of the wells, and describe them with polarization density of the form: P(r, z) = ~P n (r)5(z — z n ), where r is an in-plane 
position vector, z n represents a coordinate of the n-th well, and P„ is a surface polarization density of the respective 
well. Optical response of the system even in the presence of inhomogeneous broadening is successfully described by 
the linear dispersion theory (see, for instance, Ref. |2(]) based upon a single oscillator phenomenological expression for 
the exciton susceptibility 



where Qq and "fhom are Is exciton frequency and relaxation parameter respectively, d is a parameter describing the 
light-exciton coupling. For our purposes, however, it is more convienient to consider explicitely equations of exciton 
dynamics coupled with Maxwell equations for the electromagnetic field. In the case of waves propagating along the 
growth direction of the structure, which is considered in the present paper, only excitons with the zero in-plane wave 
number are coupled to light, and the equations of motion can be presented in the form: 
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where fi ra is the exciton frequency of the 77,-th QW. The strength of the exciton-photon coupling in Eq. (|^) is described 
by the parameter T n equal to the radiative decay rate of a single n-th well. This parameter is related to the parameter 
d of Eq. (|l|) according to d 2 = (cT)/(ttw), where w characterizes the spatial extent of exciton wave functions in the 
growth direction (we assumed it to be zero when we described the exciton polarization by the 5-function). In what 
follows, we refer to T n as to a coupling parameter. The relation of T n to the radiative life time can be established, 
for instance, if one uses Eqs. (0) and (0) to obtain the reflection coefficient for a single n-th QW, r„, 
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which in the resonance approximation (the last expression in Eq. (^)) coincides with the standard linear dispersion 
theory expression used in Refs. [2]|^ and others. 

In an infinite pure system, all r„ = Tq, Q n = f2o, z n = na, where a = Xq/2 is the Bragg's interwell separation. 
Eq. (H) describes an electromagnetic wave of one of the degenerate transverse polarizations, propagating in the 
growth direction of the structure. This equation coincides with equations used in Refs. |l6|-|l9| for one-dimensional 
chains of atoms. The spectr\i^pO|f— ideal periodic MQW's with homogeneous broadening parameter 7/j 0m = has 
been studied in many papers. E)Si3€BEiliH3 In the specific case of Bragg structures, the exciton resonance frequency, 

Oo, is at the center of the bandgap determined by the inequality uji < lo < lu u , where tui = S7o ^1 — -\/2ro/7rf2o^ 

and lo u = Qq f 1 + ^/^o/^^oj S In the formed bandgap the electromagnetic waves decay with the penetration 
(localization) length, 
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where fco — tu/c and 
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Parameter (3 is proportional to the susceptibility of the well in the absence of exciton relaxation, Eq. ([!]), and is the 
main quantity used in our theory to characterize the optical response of the single well. The homogeneous broadening 
can be taken into account by adding 2z7/ iom u> in the denominator of Eq. (^); we, however, neglect it until the last 
section of the paper, where transmission and reflection spectra of defect MQW structures are discussed. 

As it was mentioned earlier, for structures, which are longer than the penetration length, Z/ oc , for the most of the 
band-gap, consideration based upon modes of the infinite periodic structure is more appropriate. Using Eq. (|5|) 
one can show that in the case of 140 (for GaAs/AlGaAs) or 200 (for InGaAs/GaAs) wells long MQW's of the type 
considered in Refs. [|,|[ hoc > Na for 90% of the bandgap of the infinite structure characterized by the boundaries u>i 
and ojp. This bandgap is the frequency region where we look for new local states associated with the defects. In this 
paper, we consider four types of such defects. First of all, one can replace an original QW with a QW with different 
exciton frequency (fi-defect). This can be experimentally achieved by varying the composition of the semiconductor 
in the well. Another possibility is to change the coupling constant (r-defect) at one of the wells. Even though an 
experimental realization of this defect in its pure form is not straightforward, it is still methodologically important to 
consider such an idealized situation in order to be able to estimate how changes in T could effect optical properties 
of the systems with real defects. The third possibility is to perturb an interwell spacing between two wells. Here we 
distinguish two defects, which we call a- and 6-defects. The former is realized when the interwell spacing between a 
pair of wells is changed (Fig. la). It can be seen that this induces a shift in the position of all wells that follow the 
defect, making it significantly nonlocal. The 6-defect is produced when one shifts just one well keeping positions of 
the rest unchanged (Fig. lb). Experimental realizations of these two defects is simple and can be done at the sample 
growth stage. In the following section we show that each of these types affects the optical properties of the MQW 
lattice in remarkably different ways. 
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FIG. 1. Two types of interwell spacing defects. The nonlocal a-defect (a) as opposed to the local fr-defect (b). Solid bars 
represent locations of QW's in the defect lattice. The empty bars represent what would have been a perfectly ordered MQW 
lattice. 



We start with SI- and T-defects which are similar in a sense that they both introduce perturbations in the equation 
of motion, Eq. (||), which are localized at one site (the diagonal disorder). Therefore, they can be studied by the 
usual Green's function technique like the one used to treat the localized phonon states in impure crystals (see, for 
instance, Ref. [lO]). Using the polarization Green's function defined in a standard way by adding <5„ ino to Eq. (||), we 
can express the polarization of the n-th QW in terms of the polarization of the defect QW, P no as 

P n = G(n - n )P no . 

(see details of the derivation, as well as an expression for G(n — no) in Ref. |l7| ). Allowing n = no, and using expression 
for the Green function from Ref. [I?] one obtains the equation for eigenfrequencies of LPM: 
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where the function F n = (f2f — oj 2 ) / (Slf — £!q) corresponds to the fi-defect and the function F r — To/ (Ti — To) 
describes the T-defect; J7i and Ti denote respective parameters of the defect layer. Eq. (Q) is an exact consequence of 
Eqs. (|) and (§), and was first derived in Ref. for the fi-defect. It has been studied in that paper in the longwave 
approximation; it was found that the equation had one real valued solution for any f^i > Slo- We find that in the case 
of the Bragg structures, there are always two solutions for both types of defects, one below Slo and one above. This 
is a manifestation of the fact that the bandgap consists of two gaps positioned one right after the other. The above 
equations can be solved approximately when Tq <C Slo, which is the case for most materials. For the fi-defect, one 
solution demonstrates a radiative shift from the defect frequency fii, 
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while the second solution splits off the upper or lower boundary depending upon the sign of f2i — Oq : 



(2) 
def 



7T Sli 

2~ 



SI, 



Slo 



(8) 



(9) 



where one chooses lo u and "— " for Si! < ^Oj and and "+" in the opposite case. This is illustrated in Fig. 2a. It can 
be seen that the shift of from the defect exciton frequency fii is negative for Sli > Slo and positive for fii < Slo- 
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The magnitude of the shift is of the order of the coupling constant Tq, which is usually rather small, and this fact, as 

(2) 

it is shown in the next section, is crucial for the optical properties of the defect. The second local mode ojgjt lies very 
close to the edges of the bandgap, and it would be, therefore, very difficult to distinguish it from the modes making 
up the allowed bands even for negligible dissipation. 
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FIG. 2. Positions of the local modes in the bandgap (1.486 - 1.496 eV) for the S7-defect (a), a-defect (b), and 6-defect (c) as 
functions of the defect strengths. Numerical values of the exciton resonant frequency and the exciton-light coupling constant 
were taken the same as in InGaAs/GaAs structures studied in Ref. 8. 

In the case of the T-defect, one again finds two solutions of Eq. (Q): 

^def = ± 2 - r ) 2 {uj u - uji) . (io) 

which exist only for < Ti < Tq and are very close to the gap boundaries. The situation is similar to the second 
solution for the fi-defect, and one can conclude, therefore, that the T-defect would not affect significantly the optical 
spectra of the system. 

The next defect we consider, the a-defect, is shown in Fig. la. One can see that this defect differs from the other 
two in a fundamental way. An increase in an interwell distance between any two wells automatically changes the 
coordinates of an infinite number of wells: z n — na for n < rid and z„ = (b — a) + na for rid < ri, where b is the 
distance between the n-th and (n + l)-th wells. Therefore, this defect is non-local and cannot be solved by the same 
methods as the two previous cases. The best approach to this situation is to consider the system of a finite length, 
L, and match solutions for n < rid and n > rid + 1 with a solution for na < z < na + (b — a), which is schematically 
shown in Fig. 3. Eigen frequencies can be derived then considering limit L — > oo. Solutions for the finite system can 
be constructed using the transfer matrix approach. The state of the system at the m-th well is described by a two 
dimensional vector v m with components E(x m ) and (l/ko)(dE(x m )/dx): 

m 

V m = II T n V = Tv . (11) 
n=l 

The 2x2 transfer matrix r„ at the n-th well is 

~ _ ( cos(fcofln) + /3sin(fco«n) sin(fcoa n ) 
" 1 — sin(fcoa„) + (3 cos(fcoa„) cos(fcoan) 
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where a n — x n+ \ — x n and (3 is denned by Eq. (g). The eigen quasi-states for a finite system, n G (1, iV), can be found 
if one looks for non-trivial solutions when no wave is incident upon the system. This corresponds to the boundary 
conditions of the form 

v = (_ r . r ) and v N =(j t 

A resulting dispersion equation for the eigenfrequencies can be expressed in terms of the elements of the total transfer 
matrix T: 



(Tii + T 22 ) - i (Tia - T 21 ) = 0. (13) 

It is clear that for finite samples the solutions of this equation are complex, reflecting the fact that the modes are not 
stationary; they have lifetimes which tend to infinity only in the limit N — > oo. 



a-defect 




Evanescent mode Evanescent mode 



FIG. 3. Matching the solutions for half-infinite perfect MQW with cavity modes xj)± one can obtain the dispersion equation 
for (quasi)local modes (in finite systems) of such a system with the defect. 



Eq. ( |13| ) for the a-defect in the infinite MQW system, after some cumbersome but straightforward algebra, can be 
presented in the form 



cot(fc &) = — 



sin(fc a) - /3A_/2 
cos(koa) — A_ 



(14) 



where A_ = cot(fcoa) + /3/2 — yD sin(fco») is one of the eigenvalues of the transfer matrix, Eq. (|12|), and D = 

— 1 + /9 2 /4 + (3cot(koa). This equation, as well as in the cases of T- and O-defects, has two solutions - above and 
below IIq- These solutions can be approximated as 



J def - 



Q 



uj u - u>i 



(-1)! 5 * 1 ] sin«/2) 



1+ mo 2 J cos«/2) 



(15) 



= eg = n + (-l)[^]cosK/2) 

where ^ = b/a, and [...] denotes an integer part. Therefore, for Tq <§C and not very large £ , £ <C (f2o/To) 1//2 — 10 2 , 
and both solutions are almost periodic functions of b/a with the period of 1, as shown in Fig. 2b. These solutions 
oscillate between respective boundaries of the gap (uj u or uj{) and the exciton frequency Oq. At the integer values 
of £, one of the frequencies or becomes equal to Qq, and the other reaches lo u or w; depending upon the 
parity of £. When £ crosses an integer value, the solution passing through fi changes continously, while the second 
one experiences a jump toward the opposite gap boundary. The observable manifestations of the defect modes (for 
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instance, transmission resonances as described in the next section of this paper) vanish when the defect frequencies 
approach the gap boundaries. This jump, therefore, would manifest itself as a disappearance of the transmission peak 
near one of the gap boundaries, and a gradual reappearance at the opposite edge as £ changes through an integer 
value. It is interesting to note again that the exciton resonance frequency Qq, which formally lies at the center of the 
gap, behaves as one of the gap boundaries. This is another manifestation of the fact that the polariton gap in the 
Bragg structures is formed by two adjacent gaps with a degenerate boundary at re- 
calculations of the defect frequencies for the 6-defect can be done in the framework of both schemes. The transfer- 
matrix approach, however, turns out to be less cumbersome. The dispersion equation for LPM's in this case takes 
the form 
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This equation can be solved approximately assuming that the splitting of the solution from the Bragg frequency, f^o, 
is much smaller than the width of the gap. Expanding different terms of Eq. (jlT]) in terms of powers of [u — f2o) and 
keeping the lowest non-zero contribution, one obtains two different solutions for frequencies of LPM's: 
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The relative splitting of these modes from Qq is of the order of (To/flo) 3 ^ 4 , and is much smaller than the relative 
width of the gap, which is proportional to (Tq/Qq) 1 / 2 in accord with our initial assumption. Similar to the case of 
the a-defect, these frequencies change periodically with £, but unlike the previous case, they both split off the center 
of the gap, fio, and at integer values of £ merge back to fio- The maximum deviations of the local frequencies from 
Oo are of the order of (r/f2o) 3 ^ 4 , they take place for half-integer values of £. 

One can notice that the 6-defect involves two adjacent wells, and in the Green's function approach one would have 
to deal with a system of two coupled equations. Accordingly, it can be expected that there must be four possible 
local frequencies (two in each half of the gap), while we found only two of them. The reason for this is that the 
Bragg condition makes a pair of frequencies, which are both above or below £lo, nearly degenerate, and the difference 
between them is much smaller than the terms we kept in our approximate solution. This can be understood if one 
notices that the transfer matrices in Eq. ( |l2] ) contain factors cos(kb), sin(fc6) and cos[fc(2a — b)], sin[fc(2a — b)] for 
the first and the second of the involved wells, respectively. Exactly at the Bragg frequency, ka — ir, these factors 
coincide leading to the degeneracy. Since the shift of the actual local frequency from the Bragg frequency is relatively 
small, the solutions remain nearly degenerate. The fact that obtained solutions are symmetrical with respect to the 
replacement of b with 2a — b is not surprising and reflects the symmetry of the transfer matrices discussed above. 



III. DEFECT LPM'S AND TRANSMITTANCE AND REFLECTANCE EXPERIMENTS 



In this Section we study how the local defect modes obtained above affect the reflection and transmission spectra 
of the finite size periodic Bragg MQW's in the presence of homogeneous broadening. Transmission and reflection 
coefficients can be expressed in terms of the elements of the total transfer matrix defined by Eq. (JTl]) as 



T=\tr = 



2detT 



(Th + T 22 ) - i (T 12 - T 2 i) 



(19) 



R=\rr = 



- T 22 ) + i (T12 + T21) 



(Tn + T 22 ) - i (T12 - T21) 



(20) 



BOh can be shown to add up to unity. In the denominators 
3h that was obtained in the previous section by matching 



Without absorption, T and R in the form of Eqs. ( |l9| ) and 
of T and R one can recognize the dispersion relation Eq. ( 
the decaying solutions on both sides of the defect. 

V- and fi-defects differ from the original QW's only in the definition of the parameter j3, therefore, they both can 
be dealt with at the same time. The exact expression for T is cumbersome, however, if the length of the system is 
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longer than the penetration length over the frequency region of interest, it can be simplified. In this case the smaller 
eigenvalue of the total transfer matrix is proportional to exp(— N a / Ii oc (lo)) , where liap-js given by Eq. (||), and can 
be neglected. The expression for the transmission coefficient then can be presented ast3 



t = 



to 



(1 



i exp (—ikoNa) e^to cosh [(N — 2hq + 1)ko] 



(21) 



where we introduce the defect parameter e = (fidef — 0) /2v-D, which is equal to zero when f3def = P\ is the 
position of the defect QW, k = l/k oc (ui), and $ = /3/(sin(fc a)v A D). For e = 0, Eq. (|2l| ) gives the transmission 
coefficient, to, of the pure system, 



in 



2e lkL exp {-kL) 
l + i[2- 0cot{ka)]/VD 



(22) 



exhibiting an exponential decay characteristic of the evanescent modes from a band gap. 

Eq. ( pi| ) describes the resonance tunneling of the electromagnetic waves through MQW with the defect. The 
equation 

1 + e = 

can be shown to coincide with the dispersion equations for local defect modes of the infinite structure in the case 
of ti and r defects. Given this fact and the structure of Eq. (21), it seems natural to assume, as we did in our 
previous papersJiS that the transmission reaches its maximum value at the frequencies of the local modes. More 
careful consideration showed, however, that the systems under consideration behave in a less trivial way, and that 
maxima of the transmission occur at frequencies shifted from the frequency of the local modes. This is quite unusual 
behavior that distinguishes the systems under consideration from other instances of resonance tunneling. Moreover, 
we shall show that contrary to our previous result ,E3 the transmission at the maximum is always equal to unity (in 
the absence of absorption). 



Evaluation of the exact analytical expression for the transmission coefficient, Eq. (19), obtained from the total 
transfer matrix, showed that the distinctive resonance transmission occurs only if local modes lie not too close to the 
boundaries or the center (which is, strictly speaking, also a degenerate boundary) of the gap. T defects, therefore, can 
be excluded from consideration, as well as one of the solutions for the O-defect. For frequencies close to the remaining 
local mode (the solution for the ^-defect given by Eq. (||) the expression for the transmission coefficient Eq. (|2l]) can 
be simplified. Expansion in the vicinity of the local mode gives 
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(23) 



Parameters Ai, A 2 and lot arc defined by 
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A (^de/ ~ ti )y/(LO def - UJi)(lO u - LO d If) , , , s 

A 2 = t t-. r tanh UN — 2nn + Ijna \, (25) 

(LO + - LOdef ) (Udef ~ U-) 

I _ (^+ ~ U d ef)(udef - LO-) j^def - tig) —nNa t0R ^ 

lot = LOdef + p: 7= =e , (2b) 

S '0 ^ (LOdef - UJi)(LO u - LOdef ) 

where lo± = ti,Q ± (lo u — 57o) ]\[2. The resonance transmission occurs when the defect layer is located in the center 

of the system N — 2hq + 1 = 0. In this case the coefficient A 2 becomes zero and Eq. ( p3| ) can be presented in the 
following form 



4 7 p (lq-ut + Q) 2 

H2 (u;-LOT) 2 +^r 
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where Q = lot — and the parameter 70 is given by 



70 = nflo 



LOdef — ^0 



-uNa 



(28) 



The transmission spectrum described by Eq. ( |27| ) has a shape known as a Fano resonancejiij where lot is the resonance 
frequency, at which the transmission turns to unity, and parameters 7^ and Q describe the width and the asymmetry 
of the resonance respectively. One can see from Eq. (|2^) that in general the transmission resonance frequency is 
shifted with respect to the frequency of the local mode. The shift, though exponentially small for long systems 
considered here, is of the same order of magnitude as the width of the resonance 70, and is, therefore, significant. 
These two frequencies, lot and u>def, coincide in the special case when u>dej = w±. The fact that the transmission is 
equal to one, in this particular case, was obtained in our previous paper Ref. 0. 




1.49246 1.49247 1.49248 1.49249 1.49250 
co, eV 



FIG. 4. The shape of the transmission maximum for 3 lengths N = 101, 201, 301. For all lengths, the O-defect is in the 
center of the system, fii = 1.492491 eV is the exciton frequency of the defect layer. The exciton resonant frequency Qo is 
1.491 eV, the exciton-light coupling constant To is 27 fieV as in the InGaAs/GaAs structures studied in Ref. 8, however, no 
homogeneous broadening is assumed. 

At to — lot ~ Q = Q\ the transmission equals zero, which is a signature of Fano resonances. Usually the presence 
of Fano resonances is associated with interaction between a discrete level and a continuum of states. This is not 
the case in the situation under consideration. Zero of the transmission in our case is caused by the fact that the 
il-defect brings to the dielectric function of the system a new pole at f2i, which describes the interaction between 
electromagnetic waves and excitons of the defect well. The penetration length diverges at Oi, and transmission, 
therefore, vanishes. The interaction with excitons, at the same time, leads to the radiative shift of the frequency 
of the local mode, and hence, the frequency of the transmission resonance, away from Qi. The combination of 
these two factors is responsible for the Fano-like shape of the transmission resonance. The actual form of the Fano 
spectrum in an ideal system without absorption is determined by the interplay between the width parameter 70 and 
the asymmetry parameter Q. The former exponentially decreases with an increase of the length of the system, while 
the latter is length independent. However, the pre-exponential factor in 7c; is of the order of the exciton resonance 
frequency, while Q is of the order of T, i.e. significantly smaller. Therefore, in principle, there are two possible 
cases: ^S> Q for shorter systems [1 <C nNa <C log(Oo/<2)], and 7^ <C Q for longer systems. In the first case, the 
transmission spectrum has a distinctively Fano-like asymmetrical shape, while in the second limit the spectrum attains 
the symmetrical Lorentzian shape characteristic for Wigner-BreightEll resonances. Fig. 4 shows the evolution of the 
shape of the transmission resonance in the absence of absorption from Fanodike to Wigner-Breightdike behavior with 
increase of the length of the system. One can also see from this figure how the position of the transmission maximum 
moves with an increase of the length. An actual possibility to observe the Fano resonance in the considered situation 
depends strongly upon the strength of absorption in the system, which must be at least smaller than Q. More detailed 
discussion of absorption related effects is given below in the present section. 

Calculation of the transmission coefficient for the a-defect can be carried out in a similar manner. Dropping the 
exponentially small contribution from the smaller eigenvalue the transmission matrix can be written in the following 
form: 



(1 + £ ) + ie~ l = Wa $ 2 sin(^(& - a))t cosh [(N - 2n )«;a] ' 
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where 



sm(^) - Sirica) - A + (l - sin(^(6 - a)) 



sin(^a) 

Similar to the previous defects, 1 4- e a — coincides with the dispersion equation for an infinite system, but the 
transmission resonance is shifted with respect to the local mode. In the vicinity of the resonance, Eq. ( |2S| ) can be 
presented as: 



T 



1 



1 + A X 



uj — lot cosh ((N — 2no)na) (1 + i ■ Ai) 



LOdef — LOT 



(30) 



where parameters A\ and Ai are again given by Eqs. (gj) and (|25|), respectively. One only needs to replace N — 2nn + l 
in Eq. (^5|) with N — 2uq, which reflects the new symmetry of the system. The frequencies of local modes LOdef are 
now given by Eqs. ([t5]) and (|l6|), and lot is defined by 



LOT — LOdef + %(LOdef — ^o) 



(UJ + - LOdef) {LOdef ~ LO-)(uJdef ~ Loj)(uj u ~ LOdef) 
{lO u ~ «0) 4 



iNa 



(31) 



The resonance transmission again occurs when N — 2uq = 0, which requires an even number of wells in the system. 
Eq. (J30|) in this case takes the standard Wigner-Breight shape 



T = 



{LO-LOTY+^r 

with the half- width, j a , given now by 

t (LO d ef ~ tto) 2 (LO d ef ~ UO if / 2 (uj u - LOdeff' 2 



la 



iNa 



{lo u - fi ) 



(32) 



(33) 



The frequency lot, where the transmission coefficient takes the maximum value of unity, is again shifted from the 
frequency of the defect mode. The two frequencies coincide, however, when the defect frequency is made equal to 
lo±, which are the same frequencies at which transmission resonance and the local mode coincide for the f2-defect. As 
one can see from Eq. (jlj), these conditions can be satisfied simultaneously for both defect frequencies of the a-defect 
when b ~ (integer + l/2)a (see Fig. 2b). In this case the position of the transmission resonance becomes independent 
of the length of the system. 

Finally, the 6-defect gives an expression for the transmission coefficient very similar to Eq. (j2l]) with the only 
distinction that e has to be replaced by a different expression, which is too cumbersome to be displayed here. The 
maximum transmission for a given defect is again achieved when the defect is in the center of the system (N is odd) . 
Expanding the transmission coefficient near the frequency of the respective local mode, one obtains 



T = 



1 



1 

1+ 4 



LO u — ^0 



LOdef 



fir 



LOT 



LOdef — LOT 



+ i ■ sinh ((JV - 2n a + l)na) 



(34) 



with lot now given by 



LOT = LOdef - 2(LO de f ~ ^o)e 



-kNo 



(35) 



Unlike other types of defects, in this case the transmission resonance is always different from LOdef- When the defect 
is at the resonance position N — 2no + 1 = 0, the transmission spectrum again takes the Wigner-Breight shape with 
the resonance width determined by parameter 7^, where jb is given by 



76 = 2 (^f-n ) 2 e _ 

LOu — On 



uNa 



(36) 



In real systems, enhancement of the transmission coefficient is usually limited by homogeneous broadening of 
exciton resonances. Two cases are possible when exciton damping is taken into account. It can suppress the resonance 
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transmission, and the presence of the local states can only be observed in not very long systems as a small enhancement 
of absorption at the local frequency. This can be called a weak coupling regime for LPM, when incident radiation 
is resonantly absorbed by local exciton states. The opposite case, when the resonance transmission persists in the 
presence of damping, can be called a strong coupling regime. In this case, there is a coherent coupling between the 
exciton and the electromagnetic field, so that the local states can be suitably called local polaritons. Qualitatively we 
can assess the effect of absorption on resonances caused by different defects by looking at the widths of the respective 
spectra. For all types of defects the width of respective resonances exponentially decreases with the length of the 
system, consequently, in sufficiently long systems all resonances dissappear. However, pre-exponential factors make 
different defects behave differently at intermediate distances. A simple qualitative estimate would require that the 
width of the resonances be smaller than the exciton relaxation parameter. Therefore, the resonances, where the 
pre-exponential factor of the width is considerably larger than the relaxation parameter, can be observed in the 
systems of intermediate length. On the one hand, the length must be greater than the localization length of the 
respective local mode. On the other hand, it must be small enough for the width of the resonance to remain larger 
than the exciton relaxation parameter. The Fano resonance arising in the case of the fi-defect though, requires special 
consideration since its vitality is determined by the asymmetry parameter Q rather than by the width parameter 
7n- Although the latter is determined by the large pre-exponent (of the order of the exciton resonance frequency 
flo) , the former is of the order of the light-exciton coupling constant To, which is much smaller. The Fano resonance 
will likely be washed out as soon as the relaxation rate exceeds this asymmetry parameter Q. This circumstance 
will prevent observation of the Fano resonance due to £l-defect in InGaAs/GaAs MQW's, experimentally studied in 
Rcf. U which is, to the best of our knowledge, the only system, where radiative coupling was observed for systems 
as long as 100 wells. In this system, the exciton resonance frequency, r2 ? an( i the exciton- light coupling parameter, 
r , were respectively equal tq-,57o = 1.491 eV, and T = 27 [ieV ', while the exciton relaxation parameter was 
estimated as jhom = 280 /ieVQ. However, the Bragg GaAs/AlGaAs MQW system studied in Ref. |^ offers much 
more favorable conditions for the observation of this effect. According to Ref. radiative and non-radiative rates in 
GaAs / Alo.3Gao.7As system consisting of 10 wells are equal respectively to To = 67 fxeV and -fhom = 12.6 fiev. If one 
could grow sufficiently long (60 - 80 wells or more) GaAs/AlGaAs MQW's with the similar relation between To and 
"fhom, the observation of the Fano resonance induced by the f2-defect may be experimentally feasible. 
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FIG. 5. Transmission, reflection and absorption coefficients for the Q,-type defect. The defect is placed in the center of 
the MQW with 201 quantum wells. The exciton frequency of the defect well fii is chosen such that (fii — fio)/fio = 1.003. 
Numerical values of the exciton resonant frequency and the exciton- light coupling constant were taken for InGaAs/GaAs 
structures studied in Ref. 8. Solid lines show results obtained in the absence of absorption, and dashed lines were calculated 
in the presence of experimentally observedj homogeneous broadening. 
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The pre-exponent of the resonance width in the case of the a-defect is of the order of magnitude of the gap width, 
which is proportional to ^/ToTIq. It is considerably larger than jhom for the InGaAs / GaAs MQW's and we expect, 
therefore, that it can be easily observed in readily available samples of this composition. As far as the 6-defect is 
concerned, it gives rise to extremely narrow transmission peaks, which are characterized by a pre-exponential factor 
of the order of IV It, therefore, will be washed out in InGaAs /GaAs MQW's, but can be reproduced in 60 — 80 
wells long GaAs / AlGaAs systems with the parameters similar to those reported in Ref. ^|. 




FIG. 6. Transmission, reflection and absorption coefficients for the fi-type defect. The defect is placed in the center of 
the MQW with 101 quantum wells. The exciton frequency of the defect well fii is chosen such that (fii — f2o)/^o = 1.003. 
Numerical values of the exciton resonant frequency and the exciton- light coupling constant ware taken for GaAs/ AlGaAs 
structures studied in Ref. 5. All results were obtained in the presence of experimentally observedfl homogeneous broadening. 



We compliment the qualitative arguments presented above with numerical evaluation of exact expressions for the 
transmission, and reflection coefficients, Eqs. ( |l9| ) and (pcf), using parameters of the mentioned InGaAs /GaAs and 
GaAs/ AlGaAs systems. To account for homogeneous broadening quantitatively, we add an imaginary part to the 
denominator of the parameter /?: 

p _ 4r ^ 

lu 2 - Ql + 2i-f hom Lu ' 

Fig. 5 shows transmission and reflection coefficients of the Bragg MQW lattice made of 201 quantum wells with the 
il-defect in the center with parameters corresponding to the InGaAs /GaAs system. One can see that absorption 
washes out the strong asymmetric pattern of the Fano resonance, but some remains of the resonance are still quite 
prominent, and can probably be observed in high quality samples. Fig. 6 represents the spectra of 101 well long 
system with parameters corresponding to GaAs/ AlGaAs. One can see a characteristic assymetric profile of the Fano 
resonance both in transmission and reflection. There is also a remarkably strong and narrow absorption line at the 
resonance frequency with more than 8-fold growth of absorption at the resonance. This behavior should be contrasted 
with the previous figure, where absorption spectrum is rather flat with only insignificant increase at the resonance 
frequency. We confirm, therefore, our conclusion that GaAs/ AlGaAs systems with the il-defect might have a potential 
for realization of local polariton modes in the regime of strong coupling. In the case of the 6-defect, strong exciton 
absorption characteristic for the InGaAs /GaAs structures washes out any resonance features from the spectrum. The 
resonances, however, survive for a system with parameters similar to those reported in Ref. [5] for GaAs /AlGaAs. In 
this case one has very narrow symmetric transmission, reflection and absorption lines (Fig. 7) located rather close to 
the center of the gap. It is interesting to compare resonance behavior of the 6-defect and the Sl-defect in the regime of 
strong coupling. The latter results in absorption that, though increased sharply at the resonance, is still rather weak, 
at the same time the transmission reaches in this case almost 80%. In the former case, however, the transmission does 
not grow that dramatically (up to about 0.4), but absorption increases by two orders of magnitude. The difference is 
caused by the different shapes of the resonances: the Fano resonance is considerably wider and, therefore, results in 
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less dramatic increase in the maximum absorption. 
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FIG. 7. Transmission, reflection and absorption coefficients for the fo-type defect. The defect is placed in the center of the 
MQW with 101 quantum wells. The defect strength is £ = b/a = 1.5. Numerical values of the exciton resonant frequency and 
the exciton- light coupling constant were taken for GaAs / AlGaAs structures studied in Ref. 5. All results were obtained in the 
presence of experimentally observed!] homogeneous broadening. 



Fig. 8 shows the transmission and reflection of the Bragg MQW lattice of 200 quantum wells with the a-defect with 
£ = b/a = 1.5 and parameters of the InGaAs/GaAs system. This is the only defect for which the strong coupling 
regime can be realized for this material. But since this is the only system, for which the samples with large number 
of wells were experimentally grown, we shall discuss the related results in more details. 
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FIG. 8. Transmission, reflection and absorption coefficients for the a-type defect. The defect is placed in the center of the 
MQW with 200 quantum wells. The defect strength is £ = b/a — 1.5. Numerical values of the exciton resonant frequency 
and the exciton- light coupling constant were taken for InGaAs/GaAs structures studied in Ref. 8. Solid lines show results 
obtained in the absence of absorption, and dashed lines were calculated in the presence of experimentally observedB homogeneous 
broadening. 
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For the a-defect with the strength £ = 1.5, one has two symmetric with respect to the center resonances. The peaks 
are very pronounced in the tranmission, and the contrast in reflection is very large - about 30%, and can be made 
even larger if one increases the length of the system. The spectra shown in Fig. 8 can, however, be affected by the 
presence of the additional resonance due to 2s excitons. In Ref. || it was shown that this resonance causes asymmetry 
of the reflection spectrum. We expect that it could have the similar effect upon spectra presented in Fig. 8 causing 
an assymetry between low- and high- frequency portions of the spectra. Even if 2s excitons contribute additional 
resonances to the spectra, they would be rather weak,Ei and could be clearly distinguished from the effects considered 
in the present work. As an additional identification tool, one can use the dependence of the defect induced features 
upon the defect parameter b/a, while a contribution from the 2s exciton resonance would not be affected by changes 
in this parameter. Besides, changing the parameters of the defect, one can move the transmission resonances closer 
to the center of the gap and farther away from the 2s exciton resonance. 

We can conclude, therefore, that the strong coupling between local excitons and local photons can be experimen- 
tally realized in readily available samples of InGaAs/GaAs MQW's. If, however, one could grow long systems of 
GaAs/AlGaAs MQW's with parameters close to those reported in Ref. ||, the defect induced optical resonances would 
be even more pronounced with more dramatic increase in transmission and absorption at the resonance frequencies. 
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FIG. 9. The distribution of the electric field at the frequency ujt in the system with the a-defect placed exactly in the middle 
without (a) and with (b) experimentally observedQ homogenous broadening. The defect strength is £ = b/a — 1.5. Numerical 
values of the exciton resonant frequency and the exciton-light coupling constant were taken for InGaAs/GaAs structures 
studied in Ref. 8. Four curves correspond to different sizes of the system: 100, 200, 400 and 1000 QW's. 

It is interesting to note that the effects of local modes on the resulting absorption of the incident wave is very weak 
in the case of weak coupling. Even though there is some small enhancement of absorption at the resonance frequencies, 
it is much weaker than absorption peaks when strong coupling is realized. In other words, local polaritons in the 
regime of strong coupling demonstrate resonance behavior in both transmission and absorption at the same time, 
while in the weak coupling regime there is only a small effect of the local states upon all optical spectra of the system. 
The explanation for this behavior lies in the spatial distribution of the electromagnetic wave intensity throughout the 
system. In the absence of absorption, the electric field at the resonance frequency decays exponentially WftW from the 
defect layer, i.e. there occurs a strong exponential enhancement of the incident field at the defect layer .E3t3 We used 
self-embedding technique, adopted for the discrete systems in Ref. to calculate numerically the electromagnetic 
field inside MQW structure. In the absence of absorption we observed the exponential increase of the field compared 
to the amplitude of the incident electromagnetic wave Ei n 

\E max \ = \E in \-e^ Ka . (37) 

Non-radiative broadening suppresses not only the resonance transmission but also the exponential enhancement of 
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the electric field. If the strong coupling regime is not realized, the intensity of the wave decreases exponentially 
throughout the sample almost as it would in the absence of the defect, and is just slightly larger at the resonance 
frequency than off-resonance. Therefore, the peaks in absorption in these cases are also only minute. At the same 
time, Fig. 9 shows that even in the presence of absorption the a-defect in the InGaAs/GaAs MQW's, which remains 
in the regime of strong coupling, demonstrates more than a three-fold enhancement of \E\ 2 at the location of the 
defect for the system of an optimal length. Fig. 10a shows the evolution of the electric field at frequency lot at the 
location of the defect as the size of the system grows (for systems with different sizes the defect is always located in 
the center of the structure) . One can see that the exponential growth of the field in the absence of absorption (solid 
curve) changes to a nonmonotonic behavior when absorption is included (dashed line). For the particular system 
under consideration, the field reaches its maximum at about iV m =450, where we see the crossover from the resonant 
enhancement regime (N < N m ) to the exponential decay (N > N m ). Enhancement of the field at the defect explains 
both transmission and absorption resonances. 




Size 

FIG. 10. (a) The dependence of the strength of the electric field at the position of the a-defect in InGaAs/GaAs structures 
(from Fig. 9) at lot on the size of the system with (dashed line) and without (solid line) homogeneous broadening. The other 
defects did not demonstrate resonance enhancement of the electric field and are not shown here, (b) The dependence of the 
strength of the electric field at the position of the a-, b-, and fi-defects in GaAs/AlGaAs structures at lot on the size of the 
system is shown. Numerical values of the exciton resonant frequency and thpi exciton-light coupling constant were taken from 
Ref. 5. All results were obtained in the presence of experimentally observedtl homogeneous broadening. The defect strengths 
are £ = b/a — 1.5 for a- and b-defects, and (fii — f2o)/flo = 0.9963 for the fi-defect. In the presence of the realistic homogeneous 
broadening, all 3 types of the defects show the resonance enhancement of the electric field. 

In the case of GaAs/AlGaAs parameters, one can see from Fig. 10b much stronger enhancement of the field at the 
resonance for three types of defects. Judging by the strength of the absorption line shown in Fig. 7, one could expect 
an increase in the field at the 6-defect by about two orders of magnitude. From Fig. 10b we see that for N = 100 
the intensity of the electric field is, indeed, enhanced by the factor of 50. The difference in absorption between b- and 
fZ-defects also signals about the difference in the field distribution. The fact that the fi-defect leads to the vanishing 
transmission at the frequency f2i translates, as computer simulations show, into the vanishing electric field everywhere 
behind the defect at this frequency - the defect acts as an almost perfect reflector. The growth of the intensity of the 
field at spectrally close frequency f2o is also limited to only 1.5 of the incident intensity, which is much smaller than 
in the case of the ^-defect. The strong enhancement of the field at the location of the 6-defect has simple physical 
explanation. The frequency of the local states due to the 6-defect lies close to the exciton frequency that is tuned 
such that A/2 = a. This means that at the frequency of the local mode all QWs but one (the defect one, see also 
Fig. lb) lie in the nodes of the electromagnetic waved Therefore, one could expect stronger then usual grouth of the 
electric field at the defect. From this considerations, one should expect the maximum effect to occur when the defect 
QW is placed in an antinode, where the electric field is the largest. We already encountered this condition, when we 
studied the transmittance in the presence of this type of defect - the most favorable value of the defect strength was 
b/a = 0.5 or b/a — 1.5. This peculiar spatial distribution of the electric field explains the narrow transmission peak 
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and strong enhancement of the electric field on the defect. This fact has great potential for applications if this type 
of the defect is realized experimentally. 

IV. CONCLUSION 

The main objective of the present paper was to demonstrate that introducing "defects" in periodic Bragg MQW 
systems, one can obtain new opportunities to tailor optical properties of quantum heterostructures. By "defects" we 
mean wells or barriers with properties different from those of the host structure. We consider four different types 
of defects and show that each of them effects the optical spectra of the structure differently. Using available from 
literature parameters of GaAs / AlGaAsu and InGaAs/ GaAsu MQW's we generate the reflection and transmission 
spectra and show that they are significantly modified in the presence of the defect layers. 

Among four types of defect considered in the paper, three are of special interest. First of them, the fi-defect, arises 
when one of the QW is replaced by a well with a different exciton resonance frequency f2i. In this case two new 
local modes arise within the band-gap of the original structure, but only one of them, which manifests radiative shift 
from fii, can significantly effect optical spectra of the system. We found that the resonance transmission of radiation 
due to the local state results in strongly assymetric Fano-like transmission spectrum. This transmission resonance 
turned out to be very sensitive to the presence of non-radiative broadening, it can only survive if the total broadening 
(including both homogeneous and non-homogeneous contributions) does not exceed the strength of radiative coupling 
between excitons and light. Using available data we concluded that the Fano resonance could not be observed in 
InGaAs/GaAs systems, but would survive in GaAs/AlGaAs MQW's provided that one could grow 60 — 100 wells 
long systems with the parameters similar to those reported in Ref. || for 10 wells long structures. MQW's with these 
parameters would demonstrate strong assymetric transmission/reflection lines with possible increase in transmission 
up to 80%. Absorption spectrum of such structure is also characterized by a strong a narrow line with absorption 
increased at resonance, but is still remanins relatively weak. The weakness of the absoprtion correlates with the fact 
that the intensity of the field at this defect increases only slightly. 

Another defect with interesting properties considered in the paper is the 6-defect, which arises when one shifts the 
position of one of the wells keeping coordinates of all others intact. This corresponds to changes in the thicknesses 
of two adjacent barriers. This defect produces two local modes with frequencies close to the center of the gap. 
Respective transmission resonances have regular symmetric Lorentzian shapes but with rather narrow widths. The 
last circumstance makes these defects also very vulnerable to the non-radiative broadening, and they can only survive 
in systems with the radiative decay rate greater than the non-radiative one (GaAs/AlGaAs may be an example of 
such systemslj). If, however, this resonance is realized, one can obtain a narrow absorption line with almost two orders 
of magnitude increase in absorption at the resonance, or what is even more important, 50-fold increase in the intensity 
of the electromagnetic field at the defect QW. 

The most robust of the local polariton modes is produced by the et-defect. In this case the width of one of the 
barriers is increased, while all other barriers remain the same. The defect produced in this way resembles a regular 
microcavity, but with optically active mirrors. This fact distinguishes the system considered here from regular cavities 
with, for instance, Bragg distributed mirrors. Local states arising in this case give rise to rather wide resonances in 
optical spectra, which can survive strong enough non-radiative broadening and can be observed in readily available 
InGaAs/GaAs samples. 

The results obtained demonstrate that one has a great variety of opportunities to tailor optical properties using 
defect MQW's. These systems can have a range of different applications. For example, the great sensitivity of the 
defect induced features of optical spectra to characteristics of the system can be used for characterization of both host 
QW's and defect wells. This sensitivity also leads to possibilities to tune the frequencies of local modes by means 
of application of stress, electric or magnetic field, or other stimuli. Another important feature with great potential 
for applications is the increase of the field intensity in the vicinity of the defect well. It can be utilized in order to 
enhance optical non-linearity of the system. It is also worth noting that the local modes considered here can be useful 
as tunable sources of narrow luminescent lines at the frequencies of the local modes. 

Besides opportunities for applications, the effects considered in the paper are interesting from the academic point 
of view. The tunneling resonances associated with the local polaritons have a number of peculiarities compared to 
other examples of resonance tunneling phenomena. For instance, the Fano-like transmission produced by the f2-defect 
appears here under new circumstances, specific for this particular system. Another interesting feature of the local 
polaritons is that the frequencies of transmission resonances are always shifted with respect to the eigen frequencies 
of the modes and depends upon the length of the system. 

Concluding, the structures considered in the present paper demonstrate a number of interesting optical effects and 
have a potential for a variety of applications. We hope that this paper would stimulate experimental observation and 
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utilization of the predicted here effects. 
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